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ABSTRACT 

We study large-scale winds driven from uniformly bright self-gravitating disks radiating near the Eddington 
limit. We show that the ratio of the radiation pressure force to the gravitational force increases with height 
to a maximum of twice its value at the disk surface. Thus, uniformly bright self-gravitating disks radiating 
at the Eddington limit are fundamentally unstable to driving large-scale winds. These results contrast with 
the spherically symmetric case, where super-Eddington luminosities are required for wind formation. We 
apply this theory to galactic winds from starburst galaxies that approach the Eddington limit for dust. For 
hydrodynamically coupled gas and dust, we find that the asymptotic velocity of the wind is Vqo — 3 (vrot) and 
that Voo oc SFR' , where (vrot) is the mean disk rotation velocity and SFR is the star formation rate, both 
of which are in agreement with observations. However, these results of the model neglect the gravitational 
potential of the surrounding dark matter halo and a (potentially massive) old stellar bulge, which both act 
to decrease Voo. A more realistic treatment shows that the flow can either be unbound, or bound, forming 
a "fountain flow" with a typical turning timescale of ftum ^0.1-1 Gyr, depending on the ratio of the mass 
and radius of the starburst disk relative to the total mass and break (or scale) radius of the dark matter halo 
or bulge. We provide quantitative criteria and scaling relations for assessing whether or not a starburst of 
given properties can drive unbound flows via the mechanism described in this paper Importantly, we note 
that because ftum is longer than the star formation timescale (gas mass/star formation rate) in the starbursts and 
ultra-luminous infrared galaxies for which our theory is most applicable, if starbursts are selected as such, they 
may be observed to have strong outflows along the line of sight with a maximum velocity Vmax comparable to 
^ 3(viot), even though their winds are in fact bound on large scales. 

Subject headings: G&laxies: Starburst — Galaxies: Formation — Galaxies: IGM — Galaxies: Halos 



L INTRODUCTION 



Galactic-scale winds are ubiquitous in stai-burst galaxies 
in bo th the local and hi^h-red sh ift universe (Heckman et al. 
1990t iHec kman et al. 2000; Pettini et al. 2001, 2002; 
Shaplev et al...2003; .Rupke et al...2005; .Sawicki et al...2008.) . 



They are important for determining the chem ical evolution 
of galaxies and the mass-metallicity relation ('Dekel & Silkl 
1986; Tremonti et al. 2004; Erb et al. 2006; Finlator & Dave 
2008t IPeeples & Shankai1l201 1 ). and as a primary source of 
metal s in the intergalactic medium (IGM; e.g., A guirre et al.. 
1200 Ih . Moreover, galactic winds are perhaps the most ex- 
treme manifestation of the feedback between star formation 
in a galaxy and its interstellar medium (ISM). This feed- 
back mechanism is crucial for understanding ga laxy forma- 
tion and evolution over cosm ic tim e ( Springel & Hernguist 
2003; Oppenheimer & Dave J2006t lOppenheimer & Dave 
"oOStiOppenheimer et al.ii201(£ 
The most well-developed model for galactic winds 

from starbursts is flie supernova-driven model of 

IChevalier & Clega (Il985h . which assumes that the energy 
from multiple stellar winds and core-collapse supernovae 
in the starburst is efficiently thermalized. The resulting hot 
flow drives gas out of the host, swe e ping up the cool ISM 
jHeckman, Lehnert & Armus" 1993; Strickland & Stevens] 
2000; Strickland et al. 2002; Stric kland & Heckman 20091 
Fuiita et aLl l2009l) . Although this model is successful in 
explaining the X-ray properties of starbursts, the recent 
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observational results that galaxies with higher star formation 
rates (SFRs) accelerate the absorbing cold gas clouds to 
higher velocities (v ocSFR*'^^) and that the wind velocity is 
correlated with the galaxy escape velocity may challenge 
the traditional hypothesis that the cool gas is accelerated by 
the ram pressure of the hot supernova-heated wind, whose 
X-ray emission temperature varies little with SFR, circular 
velocity, and host galaxy mass, indicating a critical galaxy 
mass below which most of the hot wind escapes. These 
new observations may instead favor momentum-driven or 
radiation pres s ure-driven mo d els for the wi nd physics (e.g., 
iMartinI 119991 iMartinI l2005l ISMner et al.l '2009; Murray, 
Quataert & Thompson 2005 [hereafter MQT05]). 

The model that galactic winds may be driven by momen- 
tum deposition provided by radiation pressure from the con- 
tinuum absorption and scattering of starlight on dust grains 
was developed by MQT05. However, the conclusions of 
MQT05 are based on an assumed isothermal potential and 
spherical geometry, and are thus most appropriate for bright 
elliptical/spheroidal galaxies in formation. On the other 
hand, the theory of radiation-driven winds from accretion 
disks from the stellar to galactic scales has also been stud- 
ied (e .g., Taiima & Fukue 1996, 1998; Prog a. Stone & Drew! 
1 1 9981 fl 999; ' ,Proga,200a ,2003,) . but none of these works con- 
sidered radiation from self-gravitating disks. 

In this paper we answer the question of whether or not 
large-scale winds can be driven by radiation pressure from 
self-gravitating disks radiating near the Eddington limit. 
These considerations are motivated by the work of Thomp- 
son, Quataert & Murray (2005; TQM05), who argued that 
radiation pressure on dust is the dominant feedback mecha- 
nism in starbursts, and that in these systems star formation 
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is Eddington-hmited. For simphcity, throughout this paper 
we assume that the disk is of uniform brightness and surface 
density. In §2, we show that such disks are fundamentally un- 
stable to wind formation because the radiation pressure force 
dominates gravity in the vertical direction above the disk sur- 
face. This result is qualitatively different from the well-known 
case in spherical symmetry. In §2, we also discuss the ap- 
plicability of this model to starbursts and then calculate the 
terminal velocity of the wind along the disk pole, and its de- 
pendence on both the SFR and galaxy escape velocity. In §3, 
we assess the importance of a spherical stellar bulge and dark 
matter halo potential. In §4, we discuss the 3-dimensional 
wind structure and estimate the total wind mass loss rate. We 
discuss our findings and conclude in §5. 

2. RADIATION-DRIVEN WINDS & THE TERMINAL VELOCITY 

We consider a disk with uniform brightness and total sur- 
face density: I(r < r^^d) = I and I](r < ro) = S, where rrad 
and ro define the outer radius of the luminous and gravitat- 
ing portion of the disk, respectively. The flux-mean opacity 
to absorption and scattering of photons is k. The gravitational 
force along the polar axis above the disk is 



/grav(z) = -27rGI] 



-27rGI] 1 - 



zrdr 



(r2 + z2)3/2 



v^^o r 



(1) 



and the vertical radiation force along the pole is (e.g., 
iProga. Stone & DrewlfToollTaiima & fSu3[T998|) 

JradU) - — I , -, , _,,, - — T . 7 ■ (2) 



c 7o (r^+z^)^ c z^ + r} 



rad 



The extra factor z/\/r^ + z^ in the radiation force integral com- 
pared with the gravitational force makes the Eddington ratio 
along the pole T{z) = |/rad(z)//grav(z)| a function of height z 



r(z) = ro 



ro 



^+^#P , (3) 

z' + dd z' + dd / 



where Fq = r(z = 0) = k//(2cGS) is the Eddington ratio at the 
disk center. A disk at the Eddington limit (Fq = 1) requires 

^Edd = 2cGYi/ K,, or flux FEdd = 2tTcGYi / K. If r^ad/fD > 1 / v2 — 

0.7, then F(z) increases along the z-axis above the disk. In 
particular, for r^d — t'D, the radiation force becomes twice the 
gravitational force as z — > oo: 



Too = F(z^>oo) = 2. 



(4) 



Because F(z) increases monotonically with z, an infinitesimal 
displacement of a test particle in the vertical direction yields 
a net vertical acceleration, and the disk is thus unstable to 
wind formation. This result for disks is qualitatively different 
from the spherical case with a central point source where F = 
/rad//grav IS Constant with radius. 

If we consider the motion of a test particle in the outflow, 
the velocity along the z-axis can be written as 



47rGzjr/5 V r 



1+z- Vl+F 



(5) 



where r = ri-^d/fo, z = z/ro, and vq is initial vertical velocity. 
The first term on the right side of equation (|5]) is the "radi- 
ation potential" along the pole, while the second term is the 
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Fig. 1. — Minimum initial velocity vq/vv to drive unbound winds (thin 
lines), and the critical height z* turning the gas from deceleration to acceler- 
ations (thick lines) as a function of Eddington ratio To in the sub-Eddington 
case 2/(7rr) < Fq < 1, where f = riad/''D=l, 0.9 and 0.8. 

gravitational potential. The right side is always positive if 
Fo > 1 and rFo > 2/tt w 0.64, and thus the gas can be ac- 
celerated to infinity. On the other hand, an unbound outflow 
is still possible in the sub-Eddington case Fq < 1 if the ini- 
tial velocity vq is sufficiently large to escape the gravitational 
potential above the disk until reaching the critical height z* 
where F(z*) = 1, because the gas is decelerated from its ini- 
tial vo until it reaches z*, it will be accelerated to infinity if 
v{z*) > 0. Using this constraint we calculate the minimum vq 
required to drive an unbound wind in the sub-Eddington case, 
and the critical height z* where the wind profile acceleration 
changes sign. Figure[T]shows the results. We introduce a char- 
acteristic velocity vv = a/^ttGStd. The calculation is applied 
for v'oo = v{z — ?> oo) > Vo, or rFo > I/tt ~ 0.64. In this case 
the wind can be finally accelerated at infinity. The important 
point here is that for a disk with sub-Eddington brightness, 
and non-zero vertical velocity vq, the gas can be first decel- 
erated and then accelerated to infinity. For the typical ran- 
dom initial velocities of gas in galaxies pv^ ^ tiGT?, we have 
vo/vc ~ {h/lrofl'^ ~ 0.22(r/r£,)'/^ where the galactic thick- 
ness scale h = 0.\r has be en assumed for the second equality 
dPownes & Solomon! 1998). For this reason we expect that 
even somewhat sub-Eddington disks can drive outflows, as 
shown in Figure [T] In contrast with the spherical case for 
which v^- Vq oc F()- 1, uniformly bright self-gravitating disks 
allow the gas above the disk to be accelerated for the case with 
Fo = 1 or even for some sub-Eddington cases. In §4, we show 
the 3-dimensional trajectories of wind particles launched from 
both Eddington and sub-Eddington disks. 

We wish to apply this simple theory to starburst galax- 
ies, which may reach the Eddington limit for dust (TQM05). 
However, this application depends on the extent to which a 
disk-like collection of point sources of radiation (stars) may 
be treated as a uniformly bright disk. In the limit that the 
dusty gas of starbursts is optically-thick to the re-radiated 
FIR emission (E ^ 0.1 - 1 g cm ; kfir ^ 1 - 10 cm^ g"'; 
see TQM05), as is reached in ULIRGs, the self-gravitating 
disks around bright AGN, and some starbursts, the approxi- 
mation of a uniform brightness disk is likely valid. However, 
for S < lO'-'g cm"2 (kuv ^ lO^'cm^ g"') the disk is opti- 
cally thin to the UV radiation of massive stars and the ap- 
proximations of this paper are invalidated. In the intermediate 
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regime 10"^ <S<0.1-lg cm"^,"* the disk is optically-thick 
to UV radiation, but optically-thin to the re-radiated FIR (the 
so-called "single-scattering" limit), the application of the sim- 
ple uniformly bright disk model is tentative, as it depends on 
the distribution of the sources relative to the dusty gas, the 
inter-star spacing relative to the vertical scale of the dusty gas, 
and the grain size distribution, which affects the overall scat- 
tering albedo. We save a detailed investigation to a future 
work, but here note two important points: (1) even for low- 
E galaxies, the fraction of diffuse (potentially scattered) UV 
light may be substantial ( Thilker et al. 2005 ), and (2) if a finite 
thickness disk radiates at the single-scattering Eddington limit 
(i.e., all photons are scattered/absorbed once in the wind; e.g., 
MQT05), and if the disk is uniformly emissive and absorp- 
tive, then, for example, kuv^uv ^ Kuv-Ptot/ruv ^ ^tot/S > 
'tFiR-PpiR, where ruv ^ kuvS and Ftot, F\jy, and F^ir are the 
total, UV, and FIR fluxes, respectivly. The radiation pressure 
is thus dominated by UV emission from the "skin" of the disk 
(the Tuv ^ 1 surface), and Fuy ^ fkdd.uv — 27rcGI]/Kuv- A. 
similar argument can be made for all other wavebands that are 
"single-scattering". 

To proceed with our application to starbursts, the character- 
istic velocity v^ is written as 



. = ^/47rGT.rD = SOOkms-^Eo^V]/^ 



\,lkpc'' 



(6) 



where we take Eq = E/1 g cm'^ = I]/(48OOM0 pc'^), and 
fD,\kpc = '"d/1 kpc. Momentarily neglecting the importance 
of the surrounding dark matter halo and the potentially mas- 
sive old stellar bulge (both of which we evaluate in §3), from 
equation (|5j with Tq = 1 and r^.^^ ~ r^, the asymptotic terminal 
velocity along the pole is 



: v,^y^T/2-lc:i SSOkms"' E^^r^^ 



(7) 



' DAkpc- 

This expression for Voo can be related to the star formation 
rate (SFR) using the Schmidt law, which relates the star for- 
mation surface density and gas surface density in galactic 
disks: Esfr oc E'j^"^ (Kennicutt 1998). We approximate Eg^s = 
0.5/g_o.5E, where fg is the gas fraction. Since SFR^ EsFR7rr|,, 
we have v^o oc E'/V]/^ ex Egi^r]/^ ex SFR^-^V-'^^' or 



Voo - 400kms-'4'o.5 



■0.5 f SFR 



VSOMgyr- 



0.36 



-0.21 
'D,\kpc^ 



(8) 



which is consistent with the observation Vqo oc sFR°-'^^''°^ 
in low-redshift ULIRGs (Martin 2005) and Voo oc SFR°-^ for 
high-stellar-rnass a nd high-SFR galaxies at redshift z ^ I 
(IWeiner et al.ll2009l: but, see Fig. 17 from Chen et al. 2010). 
The observed scatter at a given SFR may be caused by differ- 
ent r/5, fg, ri-ad/ro, vq, bulge and dark matter halo mass and 
the time dependence of wind fountain with the evolution of 
the stellar population (see §3). Since we only use a simplified 
uniform disk model to derive equations d?) and dSJ, a more 
definitive comparison with the data should await a model with 
realistic distributions of surface density, opacity, and bright- 
ness. 

If we assume the galactic disk is in radial centrifugal bal- 
ance with a Keplerian velocity Viot ^ \/G7rEr, and we take the 
typical terminal velocity from equation (|7]l to compare to the 
mean rotation velocity as half of that at the radius ro, we have 



:4v/7r/2-l( 



(9) 



This estimate is again made in the absence of galactic bulge 
and dark matter halo. That the terminal velocity of the 
wind increases linearly with the galactic rotation velocity 
is also consistent with the observational results in Martin 
(2005, Fig 7). Moreover, in their cosmological simulations of 
structure formation and IGM enrichment by galactic winds, 
lOppenheimer & Pavel (12006') assumed Voo = ^o-\/To- 1 ^ 
3a, where a is the galactic velocity dispersio n (see also 
lOppenheimer & Pavel 120081: lOppenheimer et al.1 [20 1 0) with 
the factor of 3 as an assumption essentially putting in by hand. 
Our equation (|9]l shows that it can be derived consistently if 
the wind is driven by radiation pressure from a self-gravitating 
disk radiating at or near the Eddington limit. However, as we 
will show in §3, the gravitational potential of the dark matter 
halo or a spherical old stellar bulge decreases the factor "3" in 
equation (|9]l to ^ 1-3, and can even lead to bound fountain 
flows. Thus, the factor of 3 in equation (|9|l cannot be strictly 
obtained by the mechanism explored in this paper since all 
starburst galaxies have dark matter halos. Nevertheless, equa- 
tions (jTJ - (Si are essential for developing an understanding of 
the maximum envelope of values for v^o, and its dependence 
on the observed properties of starbursts. More physics of v^o 
and its z-dependence is presented in §3. 

Finally, we note that the characteristic timescale for the 
wind to reach its asymptotic velocity is 

t, ^ V/>/(4GE) ^ 3.5 X lO^E-'/'r^/J,^^, yr, (10) 

which can be compared with the timescale of a bright star- 
forming disk: f* ^ Mgas/SFR. Again employing the Schmidt 
Law, we find that 



f.^2xioVo:5sryr 



(11) 



or tc/U ^ 0.02r 
accelerated to 
formation.^ 



1/2 



fO.4 
cJg.O. 



V-0. 1 
S'^O ■ 



Therefore, the wind can be 



DAkpC gx 

' Voo before the gas supply is depleted by star 



3. BULGE AND DARK MATTER HALO 

The galactic bulge and dark matter halo around the galactic 
disk are important to the wind dynamics. If we do not con- 
sider the luminosity from the galactic bulge,'' both the bulge 
and halo only act to decrease the asymptotic wind velocity, 
and may even cause the wind to fall back to the disk as a 
"fountain flow". For the galactic bulge we employ a trun- 
cated constant density sphere, and we adopt the NEW poten- 
tial (Navarro, Frenk & White 1996) to describe the dark mat- 
ter halo distribution Pdm('") oc BT^iR + VsY^ , where r, is the 
scale radius and R = \/r^ + z^ is the distance to the halo center 
For simplicity, in this section we take r,ad ~ ro- The Edding- 
ton limit To = 1 including the dark matter halo becomes 

TTK/ 1 GMhalo / fh\ 

= 27rGE+ — = 27rGE 1 + — , (12) 

c 2r2/(cvu) V 2c';' 

where Cvir = ^vir/r?, ''vir is the virial radius, /(cyir) = ln(l + 
Cvii-)-Cvir/(l+Cvir), and 



.//. 



Mha 



alo 



Mhalo frD\ 1 



27rr,roE/(Cvir) 2Mdisk V ''^ / /(Cvir) 



(13) 



These bound.s on E depend linearly on the gas-to-du.st ratio. 



^ For an instantaneous burst of star formation, the bolometric luminosity 
decreases on a timescale ;ms ~ 4 x 10^ yr, the mai n sequence lifetime of the 
most massive stars. Comparing with equation jlOt we see that tc < fMS only 
when S/r > 3 X 10' Mq kpc-^. 

* The case of a bright spherical bulge was considered in MQT05. 
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Fig. 2. — Contours of the asymptotic kinetic energy (voo/iv) as a function of fo = fb+fh (assuming ft, = 0) and c' = rs/ro for unbound particles (/«//"), 
and the turning point logiQiztmu/ro) along the pole in the case of bound pailicles, assuming I'o = (right). The dashed line in both panels divides "bound" and 
"unbound" trajectories. 



We introduce the parameters c' = rs/ro and 



fb = 



1 /M, 



bulge 

Mdisk 



rp 



(14) 



The parameters //, and ft measure the importance of the halo 
and bulge, respectively, in determining the dynamics of the 
flow. The asymptotic velocity in terms of these parameters is 
(see eq.|5]l 



:ro 



i.AV_ 

2c' 2 



(l+fb + .fh)- 



(15) 



Note that for/},,//, -> 0, equations ( fTSl l and d?) are equivalent. 
We combine the effects of the galactic bulge and dark mat- 
ter halo using the parameter /o = /;, + fb. Assuming To = 1 , the 
condition for matter to be unbound is 



/o< 



— + 



^-1 

2 



1-^ 
4c' 



(16) 



For 1 kpc-sized disk and a typical da rk matter halo c' = 
r^/ro - 10-50 kpc/1 kpc - 10-50 (iMerritt et all l2006h . 
and an unbound wind for such a disk with Fq = 1 requires 
/o<0.6 + (vo/vc)'. 

What is a typical value for /o? If we assume a galaxy with 
negligible bulge {ft, < //,), then 



/o^0.25 



Mhaio /Mdisk \ ( ro/r. 



30 



1/30 ; V/(Cvir) 



(17) 



where we hav e taken representative values of Mhaio/A^disk ^ 
10-100(e.g.. lLeaufliaudeTai]|20TTl). c' = r ,/rn ~ 10-100. 
and /(cvh) - 2 ("for c.,w ^ 15: Maccio et"ani2008h . We find 
that for a realistic range of paramaters, /o is less than ^ 0.6 
(eq. [T6] l, which shows that matter accelerated by radiation 



pressure may be unbound in many cases. For ultra-luminous 
infrared galaxies like Arp 220, with physical sizes of or- 
der ^ lOOpc (e.g., Downes & Solomon 1998), the factor 
c' = rs/ro may be as large as ^ lO-'. However, the ratio 
of the star-forming disk mass to the total halo mass may 
be Mdisk/Mhaio ^ lO""*, so that /;, is again ^ 0.3, potentially 
implying that in very bright systems outflows can be driven 
by radiation pressure via the mechanism described here. As 
equation ( fT3] l shows, larger /o corresponds to larger ro, or 
smaller Mdisk, for a fixed dark matter halo. 

More generally, taking /o as a parameter, the terminal ve- 
locity from equation ( fTSl l is less than the value from equation 
(|7]i by a factor 



V = 



(t 



i)f' 



1/2 



(18) 



As a result, both equations ([8]l and (|9]l in §2 need to be multi- 
plied by V when the halo or bulge potential is included. For 
the purposes of an estimate, taking c' ~ 10, the relation be- 
tween Voo and (vrot) becomes 



3(1 -1.6/0) 



1/2/ 



(19) 



which shows that the dark matter halo potential well is too 
deep for particles to escape for /o > 0.6. For /o ~ 0.25 in 
equation (fTTI i. instead of the naive estimate without the dark 
matter halo in equation (|9]l, which yields Voo — 3(vrot), we ob- 
tain v'oo — 2(v,ot). Even so, equation ( fT9] l makes it clear that 
the factor o f 3 in the relation v^, c ± 3 (v'lot) assumed in the cal- 
culation bv Oppenheimer & Pavel ( I2006.) (see discussion after 
eq. |9l) cannot be stricly achieved for Eddington-limited star- 
bursts by the mechanism discussed in this paper since all sys- 
tems have dark matter halos, a nd thus non-zero /p. To the ex - 
tent that the "3" is required bv lOppenheimer & Pavel (l2006l) . 
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Fig. 3. — Outflow velocity v-/vc along the pole with c' = 10 and /o from 
top down /o = 0,0.2,0.4,0.6,0^8, 1.5 (thick lines), f ' = 4 (dashed line) and 
c' = 100 (dotted line) with /o = 1. The horizon line labeled as Voo/iv is the 
analytical solution as in equation 0. 

additional physics such as an effectively super-Eddington 
galaxy luminosity, supernova-driven hot flows, or the mech- 
anism discussed in Murray, Menard, & Thompson (2011) 
would be necessary to add to the models here. 

The quantitative scaling relations for bound and unbound 
outflows, encapsulated in the factors /o and c', are illustrated 
in Figure |2] The left panel shows contours of the asymptotic 
kinetic energy (voo/Vf)^ as a function of /o and c' = rs/ro as- 
suming v'o = for unbound outflows.^ The right panel of Fig- 
ure |2] shows contours of logjo[ztui-n/'"z5]^ the turning point of 
the wind scaled to the disk radius for bound "outflows." All 
else being equal, a disk with a more massive dark matter halo 
has larger /o, smaller Ztum, and the flow has a shorter timescale 
for reaching the turning point ftum ^ Ztum/Vf -- tcizmm/ro) (see 
eq.fTOll. The right panel of Figure |2] shows that Ztum/ro can be 
larger than 10^, implying that the flow reaches many tens of 
kpc in height above the disk before falling back towards the 
host on a timescale of ^ Gyr. 

Importantly, even if /o is large enough that the flow is 
bound, one may still observe an outgoing wind from a star- 
burst galaxy for two reasons. First, there is a region of pa- 
rameter space where the time to reach the turning point ftum is 
larger than the time for the starburst to deplete its gas supply 
f* =M^/SFR (see eq.IlIJ- Since ftum/f* ^ {ztm-n/rD)tc/t* (see 
eq.[TOll, Ztum/ro > 60r^|(^p^,/^o j Sq ' is required for ftum/f* > 
1. Second, an observable flow along the line of sight to a 
non-edge-on disk can have a maximum velocity v^ax that is 
comparable to Vqo in equation d?), even though the flow is 
bound on large scales (ztum) by the dark matter halo. Figure 
[3]shows the one-dimensional outflow velocity v'-/vc along the 
pole with various /o and c'. In the bound case, the maximum 
outflow velocities v^ax peak around 1 - lOro before decreas- 
ing to much lower values. The maximum v^ax is still com- 
parable to Voo in the absence of a halo or bulge, except for 
very large /o > 5 (see eq. \Vf\ . The velocity of the flow be- 
fore changing sign and falling back to the host galaxy is time- 
dependent and approximately reaches its maximum at a time 
~f,. (eq.Iinil. 

' Increasing vq to 0.2iv does not change the position of the contours ap- 
preciably. 



The fact that ftum > ?* and that Vmax approaches a few x (vvot) 
together imply that if galaxies are selected as bright starbursts, 
they may appear to have unbound outflows even though the 
gas is in fact bound on large scales by the dark matter halo 
potential. Indeed, Figure[3]implies that starbursts with bound 
flows may exhibit a rough correlation of the form y,nax ^ 1 - 

3(v'rot)- 

4. 3-DIMENSIONAL WINDS AND MASS LOSS RATE 

So far we have focused on the forces along the polar z- 
axis. Figure m shows 2-dimensional (2D) projections of the 
3D orbits of test particles accelerated by radiation pressure 
and gravity, as well as their constant time surfaces above a 
uniform disk, starting from an initial height 0.1 r/j, both with 
and without an NFW potential, and with no galactic bulge. 
In the absence of a halo, the upper middle panel shows that 
winds can be driven even when the disk is sub-Eddington 
(Fo = 0.88). The lower panels show that the character of the 
flow changes significantly as /o increases from 0.2, to 0.6, to 
1.0 (left to right; compare with eq.fTTIi. For the case /o = 0.2, 
the wind is clearly unbound. For /o = 0.6, particles emerging 
near the z-axis are accelerated to very large vertical distances, 
whereas particles emerging from the outer disk region fall 
back to the disk rapidly. A more massive halo with /o =1.0 
produces only a "fountain flow" in which particles fall back to 
the disk with a timescale of ^ 0.1 - 1 Gyr. Moreover, Figure 
|5] gives the 3D orbit vertical component v, evolution of the 
particles driven from different disk regions. For the face-on 
disk, the velocity v. of the wind is just the velocity along the 
line of sight. The maximum velocities are reached roughly at 
t ^ tc- Also, the maximum and terminal velocities of particles 
from the outer disk region are smaller than those from the in- 
ner disk region. As in Figure HI matter from the outer region 
of the disk can be bound by the halo's gravitational potential 
even though matter from the inner region is unbound. 

To estimate the rate of mass ejection from Eddington- 
limited disks, we first note that the flux and luminosity are 



^Edd = IttcGYs/k ~ 3 X lO'^EoK^' L© kpc" 



and 
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respectively, where ki = (k/10) cm^ g ' is the flux-mean dust 
opacity. The total mass ejection rate from the disk surface Mej 
is a local quantity measured on a disk scale height, where the 
gravity of the surrounding bulge and halo is much weaker than 
the disk's gravitational and radiation pressure forces (as long 
as /() is not too large). Therefore Mgj is only determined by 
the disk, and it can be estimated in the absence of bulge and 
halo from the integrated momentum equation. In the single- 
scattering limit, an estimate of Mej is 

MqVoc,fo=0 ~ ^Edd/C, (22) 

where Voo./o=() is the terminal velocity of the flow without a 
bulge or dark matter halo (eq. |7]i. Since Vco,/o=o is an upper 
limit to the velocity of the flow (valid as /o —> 0; see eq. \T% . 
equation ( l22b is only approximate. Nevertheless, to the ex- 
tent that Vmax is of order Voo./o=o (see Fig. [3]), equation (l22l i 
should yield an order-of-magnitude estimate of the mass loss 
rate from the disk itself. Combining equations (jT]) and ( |22] | 
we have 

2 
^2/'"rad\ ^1/2,^-1^/2 ,^ ,,,-1 ^23) 
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Fig. 4. — The particle orbits from tlie disk without (upper panels) and with (lower panels) NFW potential. Solid lines show 2-dimensional projections of the 
3-dimensional orbits. Upper panels, left to right, we take Fq = 1,0.88 and 5. The particle positions and constant time surfaces are labeled ait/tc = 2, 4, 6, 8, 10..., 
where tc = \/''d/(4GS) (eg. 1101 . Lower panels, left to right: /o = 0.2, 0.6, 1.0 with Fq = 1 and c' = 10. Constant time surfaces are marked at tjtc—'i, 10, 15, 20, 
25.... 



which is similar to observational results (e.g., Martin 2005, 
2006). Keep in mind that k\ is between the FIR limit kfir ^ 
1-10 cm^ g"' and UV limit kuv ^ lO"* cm^ g"' . The estimate 
for Mej thus varies by 2 - 3 orders of magnitude from FIR- 
thick to UV-thin disks. Using L = ec^ SFR, where e_3 = e/10"^ 
is an IMF-dependent constant, to calculate the luminosity of a 
starburst (e.g., Kennicutt 1998), the single-scattering estimate 
for the ratio of the mass ejection rate to the SFR is 



SFR 



ec 



Voo,/o=0 



'O.Se-sS; 



-1/2 -1/2 
) ''d, Ikpc 



(24) 



which implies that radiation pressure drives more matter 
from disks in low-mass galaxies (see MQT05): for example, 
Mej /SFR > 10 for Srz5 < 3.1 X IO^'Mq pc"'. However, as 
low ma ss galaxies are usii ally the most dark matter halo dom- 
inated dPersicetal .11 19961) . with lai-ger /o (see eqs.[T6]&[T7]l. 
the matter lost from low mass disks can hardly escape the halo 
potential, and will fall back on the timescale ftum (see the right 
panel of Fig. |2]l. Otherwise, some additional physical mech- 
anism beyond that described in this paper (e.g., a supernova- 
heated wind) is needed to unbind it from the halo. 

The quantity Mej in the above expressions is the rate at 
which matter is ejected from the disk in either a bound or un- 
bound flow. In the former case, as discussed in §3, the outflow 



reaches a maximum outward velocity v^ax before returning to 
the disk on a timescale ftum, which in many cases is > 4. In 
the case of an unbound outflow, Mej is a rough estimate for 
the mass loss rate from the halo as a whole on large scales. As 
emphasized in §3, the critical parameter determining whether 
or not the flow is bound or unbound is /o (eqs.[T6l&[T7]i. 

5. CONCLUSIONS & DISCUSSION 

In this paper we study the large-scale winds from uni- 
formly bright self-gravitating disks radiating near the Edding- 
ton limit. Different from the spherical case, where the Ed- 
dington ratio T = fadl fgmv is a constant with distance from the 
source, for disks F increases to a maximum of twice its value 
at the disk surface (§2). As a result, such disks radiating at 
(or even somewhat below; see Fig.[Tl) the Eddington limit are 
unstable to driving large-scale winds by radiation pressure. 

We quantify the characteristics of the resulting outflow in 
the context of Eddington-limited starbursts, motivated by the 
work of TQM05 who argued that radiation pressure on dust 
is the dominant feedback process in starburst galaxies. We 
find that the asymptotic terminal velocity along the polar di- 
rection from disks without a stellar bulge or dark matter halo 
is v'oo ^ V^ttCSto ^ 3(v,ot), where ro is the disk radius, E 
is the surface density, Viot is the disk rotation velocity (see 
eqs. I2] and |9]l, and may range from ^ 50- 1000km s"' for 
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Fig. 5. — Left panel: The time evolution of v^/vc for r/ro = 0.2 to 1.0 starting with z/ro = 0.1 and Fq = 1. The horizon line Voo/vc is the same as Figure[3] 
Right panel: the velocity v-/vc starting from (r/rD-.z/zD) = (0.2,0.1) (thick lines) and (0.8,0.1) (thin lines) with different values /o for the dark matter potential 
withc' = 10andro = 1. 



starbursts, depending on the system considered. Furthermore, 
by employing the observed Schmidt law, we find that v^o oc 
SPrO.36j.-o.2i (-ggg gq igj These results, in the absence of 
dark matter halo and bulge, are in agreement with recent ob- 
servations (e.g., Martin 2005, 2006; Weiner et al. 2009; see 
also Chen et al. 2010). The typical mass loss rate from an 
Eddington-limited disk in the single-scattering limit is given 
by equation ( |23T l and suggests these outflows may efficiently 
remove mass from the disk (eq.l24ll. 

However, both wind velocities and outflow rates can be sig- 
nificantly decreased by the presence of a spherical old stel- 
lar bulge or dark matter halo potential (§3). Deeper or more 
extended spherical potentials cause the flow to be bound on 
large scales, and to produce only "fountain flows" where par- 
ticles fall back to the disk on a typical timescale of ~ 0.1 - 1 
Gyr, depending on the parameters of the system considered 
(see Figs. |2]|5]l. The criterion for the flow to become bound 
is given in equation ( fT6l ). For typical values of the parameter 
/o (see eq. \Vl\ . we find that the winds from starbursts can be 
either bound or unbound (see Fig.|2]l. For /o ~ 0.25 in equa- 
tion ( fTTj l. the asymptotic velocity of the wind is decreased 
by a factor of ~ 0.7 from the case neglecting the dark matter 
halo completely (compare eqs. |9l and [T9]l, from ~ 3(vrot) to 
~ 2.3(vrot). An increase in/o by a factor of 2 would decrease 
the asymptotic velocity further to ~ 1.3(vrot), whereas an in- 
crease of /() by a factor of ^ 3 to ^ 0.75 would cause the flow 
to become bound. 

Importantly, even in the limit of bound fountain flows, if the 
timescale ftum for reaching the turning point Zturn is longer than 
the lifetime of the starburst f* (eq. fTTT i. one may still observe 
an outward going wind while the starburst is active and bright. 
The maximum positive velocity of the flow y^ax along the line 



of sight for non-edge-on disks is still correlated with Viot of 
the disk (see §3). These facts may complicate the inference 
from observations of winds that they are unbound from the 
surrounding large-scale dark matter halo. 

Clearly, more detailed work is required to fully assess ra- 
diation pressure on dust as the mechanism for launching cool 
gas from starburst disks. The main limitations of the work 
presented here are that (1) we consider only disks with con- 
stant brightness and surface density, (2) we do not compute 
the hydrodynamics of the flow, but instead treat the out- 
flow in the test-particle limit, and (3) we ignore other phys- 
ical effects such as the momentum and energy input from 
supernovae (Chevalier & Clegg 1985; Strickland & Heck- 
man 2009), which may act in concert with radiation pressure 
(MQT05; Mun-ay et al. 2010). The general case, with hy- 
drodynamics, realistic disk brightness, surface density, and 
dust opacity profiles will modify the picture presented here. 
Such an effort is underway (Zhang & Thompson, in prepa- 
ration). However, the basic conclusion that uniformly bright 
self-gravitating disks radiating near the Eddington limit are 
able to drive large-scale winds — particularly in the high- 
E limit in starbursts (see §2) — should not be fundamentally 
changed by more elaborate considerations. Indeed, although 
we have specialized the discussion to starbursts and dust opac- 
ity, the instability derived in §2 is of general applicability. 
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